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Abstract 

We consider the linearized instability of 2D irrotational solitary water 
waves. The maxima of energy and the travel speed of solitary waves are 
not obtained at the highest wave, which has a 120 degree angle at the crest. 
Under the assumption of non-existence of secondary bifurcation which is 
confirmed numerically, we prove linear instability of solitary waves which 
are higher than the wave of maximal energy and lower than the wave of 
maximal travel speed. It is also shown that there exist unstable solitary 
waves approaching the highest wave. The unstable waves are of large 
amplitude and therefore this type of instability can not be captured by 
the approximate models derived under small amplitude assumptions. For 
the proof, we introduce a family of nonlocal dispersion operators to relate 
the linear instability problem with the elliptic nature of solitary waves. A 
continuity argument with a moving kernel formula is used to study these 
dispersion operators to yield the instability criterion. 



1 Introduction 

Preliminaries. The water-wave problem in its simplest form concerns two- 
dimensional motion of an incompressible inviscid liquid with a free surface, 
acted on only by gravity. Suppose, for definiteness, that in the (x, y)-Cartesian 
coordinates gravity acts in the negative y-direction and that the liquid at time 
t occupies the region bounded from above by the free surface y = rj(t; x) and 
from below by the flat bottom y — —h, where h > is the water depth. In the 
fluid region {{x,y) : —h < y < r](t;x)}, the velocity field (u(t;x,y),v(t;x,y)) 
satisfies the incompressibility condition 

d x u + d y v = (1.1) 

and the Euler equation 

\d t u + ud x u + vdyU = -d x P 
1 dtv + ud x v + vdyV = —d y P — g, 
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where P(t; x, y) is the pressure and g > denotes the gravitational constant 
of acceleration. The kinematic and dynamic boundary conditions at the free 
surface {y = r](t;x)} 

v — d t r\ + ud x r] and P = P atm (1.3) 

express, respectively, that the boundary moves with the velocity of the fluid par- 
ticles at the boundary and that the pressure at the surface equals the constant 
atmospheric pressure P a tm- The impermeability condition at the flat bottom 
states that 

v = at {y = -h}. (1.4) 

In this paper we consider the irrotational case with curlw = 0, for which the 
Euler equation (|1.1[) - (|1.4[) is reduced to the Bernoulli equation 

where <f> is the vector potential such that (u, v) = V0. The local well-posedness 
of the full water wave problem was proved by Wu ( 68J) for deep water and by 
Lannes for water of finite depth f |34|). 

We consider a traveling solitary wave solution of (|l.l|l - (|1.4|l . that is, a solu- 
tion for which the velocity field, the wave profile and the pressure have space- 
time dependence (x + ct,y), where c > is the speed of wave propagation. 
With respect to a frame of reference moving with the speed c, the wave profile 
appears to be stationary and the flow is steady. It is traditional in the traveling- 
wave problem to define the relative stream function ip(x,y) and vector potential 
4>{x,y) such that: 

ip x = -v, ip y = u + c (1.5) 

and 

<fi x = u + c, 4>y = v. (1.6) 

The boundary conditions at infinity are 

(u, v) — > (0, 0) , f] (x) — > 0, as \x\ — > +oo. 

The solitary wave problem for (|l.l[) ~ (|1.4p is then reduced to an elliptic problem 
with the free boundary {y = rj(x)} (U): 

Find r)(x) and ip(x,y), in {(x,y) : — oo < x < +oo, —h < y < r)(x)}, such 
that 

A^ = in - h<y <r)(x), (1.7a) 

ip = on y = r](x), (1.7b) 

\Vtp\ 2 +2gy = c 2 on y = r](x), (1.7c) 

ip = —ch on y = —h, (l-?d) 

with 

\7ip — * (0, c) , 77 (x) — > 0, as \x\ — > +00. 
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First we give a summary of the existence theory of solitary water waves. Lavren- 
tiev ([35]) got the first proof of the existence of small solitary waves by studying 
the long wave limit. A direct construction of small solitary waves was given by 
Fridrichs and Hyers ([IS]), and their proof was readdressed by Beale ([5]) via 
the Nash-Moser method. The existence of large amplitude solitary waves was 
shown by Amick and Toland ([I]). The highest wave was also shown to exist by 
Toland ([66]), and its angle at the crest was shown to be 120 degree (Stokes's 
Conjecture in 1880) by Amick, Toland and Fraenkel ([5]). The symmetry of 
solitary waves was studied by Craig and Sternberg ([13]). Plotnikov ([55]) stud- 
ied the secondary bifurcation and showed that solitary waves are not unique for 
certain travelling speed. The particle trajectory for solitary waves was studied 
by Constantin and Escher ([21]). We list some properties of the solitary waves 
which will be used in the study of their stability. Denote the Froude number by 




and the Nekrasov parameter by 
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where q c is the (relative) speed of the flow at the crest. We note that /i is the 
bifurcation parameter used in [4]. The highest wave corresponds to /i = +oo 
since q c = O.The following properties of solitary waves are proved: 

(PI) ([?]) There exists a curve of solitary waves that are symmetric, positive 
(r) > 0) and monotonically decay on either side of the crest, with the parameter 
[i 6 (~,+oo). When /i / +oo, the solitary waves tend to the highest wave 
with the 120 degree angle at the crest. When \i \ -, the solitary waves tend 
to the small waves constructed in [5S] and [5J. Moreover, we have Vip — ► (0, c), 
T) (x) — > exponentially as \x\ — > +oo. Below, we call this solitary wave curve 
the primary branch. 

(P2) ([61], [4], [53]) Any positive and symmetric solitary wave which decays 
monotonically on either side of its crest is supercritical, that is, F > 1 or 
equivalently c > \fgh- The limit of small waves corresponds to F \ 1 ([4], 
» 

(P3) ([13]) Any supercritical solitary wave (F > 1) is symmetric, positive 
and decays monotonically on either side of its crest. Moreover, any nontrivial 
solitary wave curve connected to the primary branch must have F > 1. 

(P4) ([58]) For small amplitudes waves with /i ss ^, there is no secondary 
bifurcation on the primary branch. When the highest wave is approached, that 
is, when [i — > +oo, there are infinitely many points on the primary branch which 
are either a secondary bifurcation point or a turning point where c' (/x) = 0. 

The property (P4) is essentially what was proved in [58], though our state- 
ment above adapts the explanation in [151 p. 245]. Moreover, numerical evi- 
dences (Q2], [33]) indicate that the following assumption holds true: 

(HI) There are no secondary bifurcation points on the primary branch. 
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Under the assumption (HI), above property (P4) implies that there are 
infinitely many turning points where c' (/i) = 0. So the travel speed c does not 
always increase with the wave amplitude, and this differs greatly from KDV 
and other approximate models for which the higher waves travel faster. More 
precisely, for full solitary water waves the travel speed obtains its maximum 
before the highest wave and then it becomes highly oscillatory near the highest 
wave. This fact was first observed from numerical computations ([7], [H]), then 
confirmed by the asymptotic analysis (05], EI])- Indeed, almost all physical 
quantities (i.e. energy and momentum) do not achieve their maxima at the 
highest wave, and are highly oscillatory around the highest wave (see above 
references). This fact turns out to imply the instability of large solitary waves, 
which was first discovered from numerical computations ([64]) and is rigorously 
proved in this paper. 

Main results. Denote by fa the first turning point where c (zz) obtains 
its global maximum and, by fa the first and also the global maximum point of 
E (/z) , where 

E{p)= f ^ ~(u 2 + v 2 )dydx+ I \gtfdx. (1.8) 

is the energy of the solitary wave with the parameter /i. Numerical computations 
(BE], EH) [52]) indicate that fa > fa, and fa is the only critical point of E (fi) 
in (%,fa)- We state it as another hypothesis: 

(H2) The energy maximum is achieved on the primary branch before the 
wave of the maximal travel speed (the first turning point). 

Theorem 1 Under the assumptions (HI) and (H2), the solitary wave is linearly 
unstable when /i G (fa, fa), where fa and fa are the maxima points of the 
travel speed and energy, respectively. The linearly instability is in the sense that 
there exists a growing mode solution e A * [q (x) , ip (x, y)] (X > 0) to the linearized 
problem A2.2$ , where r) (x) , ip (x, y) G C°° H H k for any k > 0. 

Our next theorem shows that there exist unstable solitary waves approaching 
the highest wave. 

Theorem 2 Under the assumption (HI), there exists infinitely many intervals 
Ii C (fa, +oo), (i = 1, 2, • ■ • ) with linin^oo max{/z| /i G /„} = +oo ; such that 
solitary waves with the parameter /i G Ii are linearly unstable in the sense of 
Theorem [71 

Theorem [5] suggests that the highest wave (n = oo) constructed in [B] is 
unstable. This contrasts with the stability of peaked solitary waves in some 
shallow water wave models ([22], [IB]). Numerical evidences ([M], [52]) suggest 
that solitary waves are spectrally stable when /i£ (^,fa), and linearly unstable 
when [i > fa , at least before the first few turning points where the computations 
are reliable. We note that the amplitude of the maximal energy wave with the 
parameter fa is already close to the maximal height ( 64 ) . So the unstable waves 
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proved in Theorems [T] and [3] are of large amplitude, and therefore this type of 
instability can not appear in approximate models which are derived under the 
small amplitude assumptions, such as KDV equation or Boussinesq systems. 
Numerical evidences ([55]) also suggest that this large amplitude instability can 
lead to wave breaking. Such wave breaking induced by large unstable waves had 
also been used to explain the breaking of waves approaching beaches ([IS], [56] . 
[57]). More discussions of these issues are found in Remarks [1] and [2] (Section 
5). 

The proof of Theorems [1] and [5] also has some implications for the spectral 
stability of solitary waves with fj, < fa. We note that the traveling waves of 
full water waves even with vorticity are shown (|12j. |31j ) to be always highly 
indefinite energy saddles under the constraints of constant momentum, mass 
etc. Therefore, their stability cannot be studied by the traditional method of 
proving (constrained) energy minimizers as in many model equations such as 
the KDV type equations ([9], [13]). So far there are few effective methods for 
proving nonlinear stability of energy saddles. So naturally, the first step is 
to study their spectral stability, namely, to show that there does not exist an 
exponentially growing solution to the linearized problem. The following theorem 
might be useful for this purpose. 

Theorem 3 Assume the hypothesis (HI). Suppose that there is a sequence of 
purely growing modes e A "* [rj n (x) ,ip n (x,y)] (A„ > 0) to the linearized problem 
for solitary waves with parameters {(i n }, with X n — > 0+ and \x n — > /io where 
is not a turning point, then we must have ^ (/io) = 0. 

By the above theorem, if an oscillatory instability can be excluded, that 
is, any growing mode is shown to be purely growing, then the transition of 
instability can only happen at the energy extrema or turning points. Numerical 
results in [BJ, [32] justify that the growing modes found are indeed purely 
growing for solitary waves before the first few turning points. If additionally 
the spectral stability of small solitary waves can be proved, then it follows that 
the solitary waves are spectrally stable up to the wave of maximal energy. 

Comments and ideas of the proof. First, we comments on related 
results in the literature. In W|. Saffman considered the spectral stability of 
periodic waves in deep water (Stokes waves), under perturbations of the same 
period (so called superharmonic perturbations). The picture of superharmonic 
instability of Stokes waves ([62]) is similar to that of the instability of solitary 
waves. The approach of [59j is to take the finite mode truncation of the linearized 
Hamiltonian formulation of Zakharov ([69 ) and study the eigenvalue problem 
for the matrix obtained. By assuming the existence of a sequence of purely 
growing modes with the growth rate A n — > 0+ and parameters /i„ — * /io, the 
solvability conditions are checked to the second order to show that /io must be 
an energy extremum. That is, an analogue of Theorem [3] was established in [59] 
for Stokes waves. However, the analysis in [59 is at a rather formal level. First, 
Zakharov's Hamiltonian formulation has a highly indefinite quadratic form that 
is unbounded from both below and above. This is due to the indefiniteness of the 
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energy functional of the pure gravity water wave problem ([12]) as mentioned 
before. So it is unclear how to pass the finite truncation results in 59 to the 
original water wave problem. Secondly, an implicit assumption in [59] is that the 
truncated matrix has the zero eigenvalue of geometric multiplicity 1 . It is unclear 
how to check and relate this assumption to the properties of steady waves. For 
solitary waves, the truncation approach of [59] seems difficult to apply because 
of the unbounded domain. Recently, in [32], [33], Kataoka recovered Saffman's 
formal result (or analogues of Theorem [3]) for periodic waters in water of finite 
depth and for interfacial solitary waves in a different way. The analysis of [32] , 
[3"5] is again formal and of similar nature as 59J. That is, by assuming the 
existence of purely growing modes with vanishing growth rates, the first two 
solvability conditions were checked to show that the limiting parameter is an 
energy extremum. We note that in the above papers of Kataoka and Saffman, 
the existence of a sequence of purely growing modes with vanishing growth rates 
was only assumed but never proved. Moreover, their analysis are perturbative, 
only for travelling waves near energy extrema. In this paper, we rigorously prove 
the linear instability of large solitary waves and our method is non-perturbative, 
which can apply to solitary waves far from the energy extrema. 

Below, we briefly discuss main ideas in the proof of Theorems [1] and [21 To 
avoid the issue of indefiniteness of energy functional, we do not adapt Zakharov's 
Hamiltonian formulation in terms of the vector potential 4> on the free surface 
and the wave profile r\. We use the linearized system derived in |31j . in terms 
of the infinitesimal perturbations of the wave profile r\ and the stream function 
ip restricted on the steady surface iS e .Then we further reduce this system to 
get a family of operator equations AQ (ip\s,,) = 0, where A > is the unstable 
eigenvalue to be found. The operator A* is the sum of the Dirichlct-Neumann 
operator and a bounded but nonlocal operator. The idea of above reduction is 
to relate the eigenvalue problems to the elliptic type problems for steady waves. 
The hodograph transformation is then used to get equivalent operators A x de- 
fined on the whole line. The existence of a purely growing mode is equivalent 
to find some A > such that the operator A x has a nontrivial kernel. This is 
achieved by using a continuity argument to exploit the difference of the spectra 
of A x near infinity and zero. 

The idea of introducing nonlocal dispersion operators with a continuity ar- 
gument to get instability criteria originates from our previous works ([40], [39], 
[38] ) on 2D ideal fluid and ID electrostatic plasma, which have also been ex- 
tended to galaxy dynamics [28] and 3D electromagnetic plasmas [41] , [42] . The 
new issue in the current case is the influence of the symmetry of the problem. 
More specifically, we need to understand the movement of the kernel of .4° that 
is due to the translation symmetry, under the perturbation of .4° to A x for small 
A. This is obtained in a moving kernel formula (Lemma [53}- The convergence 
of A x to .4° is very weak, so the usual perturbation theories do not apply and 
the asymptotic perturbation theory by Vock and Hunziker ([ST]) has to be used 
to study perturbations of the eigenvalues of .4°. An important technical part 
in our proof is to use the supercritical property F > 1 and the decay of solitary 
waves to obtain a priori estimates and gain certain compactness. In particular, 
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F > 1 implies that the essential spectra of the operators A x lie in the right 
half complex plane. The techniques developed in this paper have been recently 
extended to show instability of large Stokes waves ([13]) and get instability cri- 
teria for periodic and solitary waves of rather general dispersive wave equations 

(S3], 05]). 

In Lemma 14.11 we prove that the zero-limiting operator A is exactly the 
same operator used in [58] for studying the secondary bifurcation of solitary 
waves. This link is interesting and a little unexpected since our derivation of 
the operator A is totally unrelated to the formulation used in [55]. We note 
that the bifurcation results in [58 have no implications for instability of solitary 
waves. Indeed, for water wave problems, there seems to be no definite relations 
between the stability and bifurcation of travelling waves. For example, it was 
shown in |31[ Remark 4.13] that the bifurcation of nontrivial traveling water 
waves are unrelated to the exchange of stability of trivial flows. From numerical 
works ([64], [E], [33]), the exchange of instability at energy extrema for solitary 
waves does not imply any secondary bifurcation there. 

This paper is organized as follows. In Section 2, we give the formulation of 
the linearized problem and derive the nonlocal dispersion operators A x . Section 
3 is devoted to study properties of the operators A x , in particular, their essential 
spectrum. In Section 4, we apply the asymptotic perturbation theory to study 
the eigenvalues of A x for A near 0. In Section 5, we derive a moving kernel 
formula and prove the main theorems. Some important formulae are proved in 
Appendix. 

2 Formulation for linear instability 

In this Section, a solitary wave solution of (|1.7[) is held fixed, as such it serves 
as the undisturbed state about which the system (| 1 . ![) — (| 1 .4|) is linearized. The 
derivation is performed in the moving frame of references, in which the wave 
profile appears to be stationary and the flow is steady. Let us denote the 
undisturbed wave profile and relative stream function by r] e (x) and ip e (x,y), 
respectively, which satisfy the system (|1.7[> . The steady relative velocity field is 

(u e (x,y) +c,v e (x,y)) = (ip ey (x,y),-ilj ex (x,y)), 

and the steady pressure P e (x,y) is determined through 

P e (x,y) = l -c 2 - i|V^(x,y)| 2 - gy. (2.1) 

Let 

V e = {(x, y) : — oo < x < +oo , — h < y < rj e (x)} 

and 

S e = {(x,T] e (x)) : OO < X < +Oo} 

denote, respectively, the undisturbed fluid domain and the steady wave profile. 



7 



Let us denote 

(r)(t; x),u(t; x, y),v(t; x, y), P(t; x, y)) 

to be the infinitesimal perturbations of the wave profile, the velocity field and 
the pressure respectively. The stream function perturbation is ip (t;x,y), such 
that (u,v) = (ip y , —ip x ) . 

The linearized water-wave problem was derived in [31] . and it takes the 
following form in the irrotational case: 

A^ = in T> e , (2.2a) 

d t r] + d T (ip ey r)) + d T ip = on S e ; (2.2b) 

P + P eyV = on S e ; (2.2c) 

d t d n ^ + d T (^ ey d n ^)+d T P = on S e ; (2.2d) 

9^ = on {y = -h}, (2.2e) 

where 

d T f = d x f + r]e X dyf and d n f = d y f - r] ex d x f 

denote the tangential and normal derivatives of a function / (x,y) on the curve 
{y = r) e (x)}. Alternatively d T f(x) = f (x , r] e (x)) . Note that the above 
linearized system may be viewed as one for tj){t\x,y) and rj(t;x). Indeed, 
P(t; x, T] e (x)) is determined through (|2.2c[) in terms of r](t; x) and other physical 
quantities are similarly determined in terms of ip(t; x, y) and rj(t; x). 

A growing mode refers to a solution to the linearized water-wave problem 
(|2~2a|) - ([2~2e|) of the form 

{r)(t;x),iP(t;x,y)) = (e xt V (x), e xt ^(x, y)) 

and P{t\x,rj e {x)) = e xt P(x,rj e (x)) with Re A > 0. For a growing mode, the 



linearized system (I2.2j) becomes 

AV> = in V e (2.3) 
and the following boundary conditions on S e , 

^V( x ) + (^(^^ Ve(x))v( x )) = "^^(^ ^(x)), (2.4) 
P(x, T) e {x)) + P ey (x, f]e(x))n(x) = 0, (2.5) 

\lpn{x) + — (lp ey {x,Tle{x))lp n (x)) = -—P(x,T] e (x)). (2.6) 

ax ax 
We impose the following boundary condition on the flat bottom 

il>(x,-h) = 0, (2.7) 



8 



from which (|2.2cp follows. In summary, the growing-mode problem for a solitary 
water-wave is to find a nontrivial solution of (|2.3[) - (|2.7[) with Re A > 0. Below, 
we look for purely growing modes with A > and reduce the system (I2.3[l - (|2.7p 
to one single equation for ip\s e - For simplicity, here and in the sequel we identify 
ip ey (x) with ip ey (x,T] e (x)) and <p(x) with 4>(x,r] e (x)), etc. First, we introduce 
the following operator 

cX= ( A + ^ ( ^ (x)0 ) (2 ' 8) 

Note that ip ey > in T> e by the maximum principle and Hopf 's principle ( 23 ) , 
and the fact that ip ey = u e + c —> c as \x\ — >oo. Thus 

CO < ^ey(x,T) e (x)) = U e + C < C x , (2.9) 

for some constant Co, c\ > 0. Defining the operator 

~jZ i^ey{x)(t){x)) , 

we can write C x as 

C x = =M (2 10) 

\ + ViP ey {x) \ \ + VjiP ey (x)' v ' ; 

Denote -L^(iS e ) to be the ipey — weighted L 2 space on S e . Because of the bound 
(|2.9p . L^^Se) and L 2 (S e ) are norm equivalent. Note that the operator V is 
anti-symmetric on LA (S e )- 



Lemma 2.1 For A > 0, define the operator £ X ' ± : L\ (S e ) — > I/B, (iS £ ) fry 



and 



A±X> 

TTien, 

(a) The operator £ X '^ is continuous in X, 

ll fA,± L| ev (S e )^L| ea (S e ) ^ ( 2 - H ) 

When A — > 0+, £ A:± converges to strongly in L^(iS e )- 
(cj When A — > +oo, £ A,± converges to 1 strongly in L?^ (S e ). 

Proof. Denote {M a ; a £ R 1 } to be the spectral measure of the self-adjoint 

2 

V>e 

2 



operator 7?. = — £2> on L^(<S e ). Then 

2 ^- / Jll 71/f J.II2 IIWJI 2 



I^IlL 



/ 


A 




A ± ia 



d\\M a (t>\\l2^ < I d\\M aYULl M ,, NiJ 
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and (|2.1ip follows. Similarly, we get the estimate (|2.12[) . To prove (b), we take 
any £ (<S e ) and denote the function £(a) to be such that £(a) = for 
a^O and £(0) = 1. Then by the dominant convergence theorem, when A — > 0+, 



A 



d\\M a <p\ 



2 



A ± ia 

aa)d\\M a (b\\ 2 L2 = ||M {O} 0" 2 

lb--- 



Note that M{ } is the projector of to kerl?={0}. So M{ } = and 

£ X ' ± tp — > in . The proof of (c) is similar to that of (b) and we skip it. I 
By the above lemma and the bound (|2.9|) on ip ey , we have 



Lemma 2.2 For A > 0, the operator C x : L 2 {S e ) -> £ 2 (<S e ) denned 6y i fO|) 
has the following properties: 
(a) 

\\r x \\ < n 

for some constant C independent of A. 

(b) When A — > 0+, C A converges to ^ 1 ^ strongly in L 2 (S e ). 

(c) When A — ► +oo, C A converges to strongly in L 2 (S e ). 

By using the operator C A , the growing mode system (|2.3|) - ()2.7j) is reduced to 

ip n (x) +C x P ey (x)C x iJj = 0, on S e , (2.13) 

Aip = in X> e , 
ip(x,-h) = 0. 

We define the following Dirichlet- Neumann operator M e : H 1 (S e ) — > L 2 (S e ) by 

W e / = <9„^/ = (flyV/ - Vexd x 1pf) (x,7] e (x)) , 

where ?/>/ is the unique solution of the following Dirichlet problem for / £ 

H 1 (S e ) 

Aipf = in T> e , 
ipf(x, —h) = 0. 

Then the existence of a purely growing mode is reduced to find some A > such 
that the operator A x defined by 

A x e =N e +C x P ey (x)C x : H 1 (S e ) -» L 2 (S e ) (2.14) 

has a nontrivial kernel. Note that if we denote by <f>f the holomorphic conjugate 
of -0/ in 2? e , then d n ibf = 7r<Pf- This motivates us to define an analogue of the 
Hilbert transformation as in [58j . by 

(CJ) (X) = [ Af e f dx. 
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Then the operator Af e can be written as W e = -j^C e . From the definition, 
C e .f + if and / — iC e f are the boundary values on S e of some analytic functions 
in T> e . Below, we further reduce the operator A\ to one defined on the real 
line. First, we define the holomorphic mapping F : T> e — > Rx (-c/i,0) by 
F (x, y) = {(j) e (x, y) , tp e (x, y)). We denote 

(f,0 = - (<l>e,i>e) eD a = Rx (-h,0) 
c 

and define the mapping G : D — > V e by 

GK,o = (i(e.o,»u,0) = -F , - 1 (ce,o?)- 

The flat Dirichlet-Neumann operator M: H 1 (R) — » £ 2 (R) is defined by 

A/"/ = <9 ? V/|{ ? =o}, 

where ^/ is the solution of the following Dirichlet problem for / £ H 1 (R) 

Aipf = in A,, 
^/|{ f =o} = /, 

^/|{ C =-ft} = 0. 

Similarly, we define the operator C by Cf = Aff dx. Then N — and 
Cf + if or / — iCf are the boundary values on = 0} of analytic functions in 
D Q . Moreover, N is a Fourier multiplier operator with the symbol ([25]) 

n(fc) = A— . (2.15) 

v ; tanh(fc/i) v ; 

To separate the uniform flow (c, 0) , we rewrite 
Denote 

f(0 =2/1 (CO). (2.16) 

Then we can set x\ (£, 0) = Cu; by adding a proper constant to the vector 
potential cf> e . The mapping G restricted on {s = 0} induces a mapping _B : 
R — >S e defined by B (£) = (£ + Ciu, w). Denote z — x + iy and p = £ + then 

<i (0 e + «-0 e ) dp 1 c 

u e + c — iv e = = c— = c-r- = . 

dz dz 2£ 1 + a^xi + lo^yi 

So on {5 = 0}, we get 

1 + Afw w' fn , _ 

u e + c = c t — , v e = c o (2.17) 

\w\ 2 \w\ 2 v ; 
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where 

W = (1 + %ei +id 6 y x ) | {<;=0} = l + Cw' + iw' = 1 + Mw + iw' (2.18) 
and ' denotes the £— derivative. From (|2.17[) , 

(lie + c) C 



l+Afw = 



(u e + c) 2 + v 2 e ' 
and thus by (|2.9p , there exists C2 , C3 > such that 

c 2 < 1 +7Vw < c 3 . (2.19) 
We define the operator B : L 2 (S e ) — > L 2 (R) by 

(B/) (0 = / (B (0) =f(Z + Cw(0,w (0) , for any / e L 2 (5 e ). 
Since BC e = CB and = (1 + A/to) B-£, we have 

BN e B~ x = B^-CeB- 1 = —±—^-C =—L—N, 
dx 1+A/W£ {1+ATw) 



and 



BA X B^ = BNeB- 1 + (BC X B-^ BP ev B~ x (bC x B~ x \ 



1 N + C x P ey (S,)C > 



Here, 



1+N 



i + AfwdC J V 1 +^ £ ^ 

and we use i>ey(0> P^y (0 to denote ^ey (B (£)) , P ey (B (£)) etc. For A > 0, we 
define the operator A x : H 1 (R) -> L 2 (R) by 

= A/" + (1 + AAi;) C x P ey (£) C x . 

Then the existence of a purely growing mode is equivalent to find some A > 
such that the operator A x has a nontrivial kernel. 

3 Properties of the operator A x 

In this section, we study the spectral properties of the operator A x . First, we 
have the following estimate for the Dirichlct-Neumann operator A/". 
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Lemma 3.1 There exists Cq > 0, such that for any Se (0,1) and f £ (R) . 
we have 

|2 , /-» e II jtm2 



1 



W,/)>(i-«) r ||/||k + <W| 



ff2 



where 



\\f\\U = I + 



and / (fc) is the Fourier transformation of f . 
Proof. By the definition (|2~T5l) . 

fc 



/(*) 



tanh (fc/i) 
It is easy to check that the function 

h(x) 

satishes 

h(x) 



dk 



tanh(|fc| h) 



f(k) 



dk. 



tanh (xh) 



x>0 



1 h(x)-l 
- > and lim — — ^ = 1. 



So there exists K > 0, such that h (x) — j- > ^x, when x > K. Thus 



(AA/,/)> - / /(*) 



2 1 
dk+- 



\k\ f(k) 



k\>K 
2 



>(l-8)±\\ff L2 +min 



6 

6 S 1 



dk 



k\<K 



/(*) 



dk 



(1 + 1*1) /(*) 



|fe|>K 



/(*) 



(/A: 



2/i,' 2i<7i' 2 

This proves the Lemma with Co = min { ^ , , 5 } • B 

We have the following properties for the operator C x . 

Lemma 3.2 For A > 0, the operator C x : L 2 (R) — > i 2 (R) defined by 
satisfies: 
(a) 



<c, 



L 2 (R)^L 2 (R) 

/or some constant C independent of A. 

(7>J When A — > 0+, C A converges to ^ strongly in L 2 (R). 

(cj When A — > +00, C A converges to strongly in L 2 (R). 

Proof. By (|2.19p . the operator 2? and 23 _1 are bounded. Since C A 
£>C A £> _1 , the above lemma follows directly from Lemma[ 
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To simply notations, we denote &(£) = 1 + Afw and define the operators 
V 



' d ^ ey (0-) and £ A <±0(x) A 



A±X> 



The operator I? is anti-symmetric in the bipey— weighted space (R)- Sim- 

ilar to the proof of Lemma 12. 1[ we have 



Lemma 3.3 (a) For any A > 0, 



< 1, 



and 



l-E 



A,± 



< 1. 



(3.1) 
(3.2) 



(b) When A — > 0+, £ ' converges to strongly in L 



He 



(c) When A — > +oo, £ > converges to 1 strongly in 



The operator C A can be written as 



1 



A 



1 



Proposition 1 For any A > 0, we /iave 



1-f 



'-(* J )c{.|H.A*i(i-£)} 



We note that 



(3.3) 



(3.4) 



by Property (P2), so the above Proposition shows that the essential spectrum 
of A x lies on the right half plane and is away from the imaginary axis. To prove 
Proposition [1] we need the following lemmas. 



Lemma 3.4 For any u 6 H * (R), we have 
(i) For any A > 0, 



A,±„ 



i < C\\u\\ i 



(3.5) 



for some constant C independent of X. 

Below, let F (£) be a fixed bounded function that decays at infinity. Then 
(ii) Given A > 0, for any e > 0, there exists a constant C £ such that 



F£ x ^u 



<e\\u\\ i 

L 2 — 1 "fl! 



c E 



L 2 
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(Hi) For any e > 0, there exists X e > 0, such that when < A < A £ , 



<e\\u\\ i 
L 2 — 11 "h? 



(iv) For any e > 0, there exists A £ > 0, such that when A > A s , 



Fli-e 



x.± 



Proof. Proof of (i): Denote |M a ;a <E R 1 } t° ^ e the spectral measure of 
the self-adjoint operator R = —iD on LfL . For s > 0, we define the space 



H s = i u G L 



6l/>e 



R 



with the norm 



HI*. = IMIl^ + 

oipey 



R 



\\ u hi, 



\a\ 2s d\\M a u\\ 2 L2 



R 



where i? is the positive self-adjoint operator defined by J|a| s dM a . We claim 
that the norm \\-\\jj s is equivalent to the norm ||-|| ffs , for < s < 1. When 
s = 0, = L^ e and = L 2 . Since 6 and ^ey are bounded with positive 



lower bounds, 



and 



L' 2 



are equivalent. When s = 1, we have 



1 d 
6d£ 



{ipeyU) 



bipeydx 



which is clearly equivalent to ||w||^i , again due to the bounds of b and ip e y When 
< s < 1, the spaces H s (H s ) are the interpolation spaces of H°(H°) and H 1 
(if 1 ). So by the general interpolation theory (QT]), we get the equivalence of 
the norms ||'||^j s and Thus, there exists G\^C% > 0, such that 



<?iNy <\\u\\ Hi <c 2 \\u\\^ 



(3.6) 



Since R and £ > are commutable, we have 



£ x ^u 



ff2 



R 



i? 



The estimate (13.51) follows from above and (13.61). 
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Proof of (ii): Suppose otherwise, then there exists Eq > and a sequence 
{u n } £ H 5 (R) such that 



FS^i 



We normalize u n by setting 



> £n II u n II i + n 



F£ x ^u r 
1 



L 2 



i + e 2 

1. Then 

1 



L 1 



i+e 



< 



So there exists £ H 2 , such that u n — > Mm weakly in _ff 2 . Since ~^ 

strongly in L 2 , we have = 0. Thus v n = £ A,± it n converges to weakly in 
L\ By (i), 

C 

IKILi < c||u„|i 1 < — . 



'H 2 



ff 2 



Let G be a cut-off function for {|£| < R}- We write 
-F , = ^Xfl + ^(l-Xfl)=-F , i+^2. 

Then 

||*3«»|Il» < C max \F (01 KIL* < C- max |F (£)| < J, 

\£\>R Eq \Z\>R Z 

when R is chosen to be big enough. Since F\ has a compact support and 
ij2 L 2 is locally compact, so F\V n — ► strongly in L 2 . Thus, when n is 
large enough, 

3 

\\Fv n \\ L2 < \\F lVn \\ L2 + \\F 2 v n \\ L2 < 



4 

F£ x ^u r 



1. 



L- 



This is a contradiction to the fact that ||Fu„|| L2 

Proof of (iii): Suppose otherwise, then there exists £0 > and a sequence 
{u n } £ i?i (R) , A„ -> 0+, such that 



L2 >eoKiy 



Normalize u n by 



^ = 1. Then ||it n ||^i < — . Let u n — * weakly 



in 7J 2 . Then for any v £ L , we have 



0. 



because by Lemma T3.31 btp ey £ Xn ' T J — ► strongly in i 2 when A„ — * 0+. 

So £ A " !± u„ — > weakly in L 2 , and this leads to a contradiction as in the proof 
of (ii). 

Proof of (iv) is the same as that of (iii), except that we use the strong 
convergence 1 — — > when A„ — ► 00. I 
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Lemma 3.5 Consider any sequence 

{u n }eH^(R), ||ti„|| 2 = l, supp u n C {C| |f | > n} . 
Then for any complex number z with Rez < ^Sq, we have 

Re ((.4 A - z) u n ,u n ) > ^S , 



when n is large enough. Here, Sq is defined by {3.$ - 
Proof. We have 

Re ((.4 A - z) u n ,u r ) = {Afu n , u n ) - Re z + Re (bC x P ey (£) C x u n , u„) . (3.7) 
For < S < 1 (to be fixed later), by Lemma |3~T1 



{Nu n ,u n ) > (l-S)- + C 6\\u n f H i . (3.8) 



Note that by (J2H]) 



Pey\s s (l V ^e| 2 ) = ~9 - {^ex^exy + Ay^eyy) 

= ~9 + ^ey {ipexx + Vexlpexy) = -.9 + ' l Pey~T~ (V'ex) ■ 

Denote 

i , e 1/ (0 = -s + ^(05(0. a(0 = ^(^«)(0- (3-9) 

Then a(£) decays exponentially when |£| — * oo. We have 

= -g(b(l- £ A <+) — f 1 - £ A -+) — tin, U 



• i — - I, ~-n i 

6(1- £ X - + ) ail- —u n , u n 



= T 1 +T 2 . 



Denote 



6(0=6-1, c(0 = -/---• (3.10) 

Then 6, c tends to zero exponentially when |£| — > oo. The first term can be 
written as 

Ti = -g (b (l - £ A <+) -L^, (l - £ A -) J-u, 

6^(1 -U, 



g [ Hey 



— i-f A '+ 

/ > 



-—U n , (l - £ A ' ) - U r 

Wey V ' Wey 



2 
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where in the above we use the fact that the operator T> is anti-symmetric in the 
space L\^ e . In the rest of this paper, we use C to denote a generic constant in 
the estimates. By Lemma I3~!?l and the assumption that supp u n C {£| |£| > n}, 
we have 



\Tl\<9 
<9 



l-£ x 
1 

-Mr. 



1 



< y- 



bip e 



v / ip e „ 



Li„ 



Li, 



1 



ey- 



Li, 

1 



Un\ d£) / 6- \u n \ d£ 



ly 



fl (i + Cm«( 6(05(0 + 6(0 +|g(0l 
• fi + Cmax( 6(05(0 + 6(0 +|£(0I 

V C £ \>n V 



Since 



we have 



\Ti\ <g 

< c 



+ o 



blp e 



l-S 



5, £ 



i x > + c, 



c£ 



1 -£ 



4>e 



V'e 



L 2 



L 2 



1 

c- — u 7 

Vey 



L 2 



L 2 



Since 5(0 decays at infinity, by Lemma [3.41 (ii), for e > (to be fixed later), 
there exists C £ such that 



So 



c£ X ' + u r 



c£ x + 



L2 <e\\u n \\ H , + C e 



< £ IWnW i + 



a 



L 2 



4>e 



L 2 



1 

< - 

C 



c£ x '+u r 



?£ x '+ (cu n ) 



< - K 

c 



HI 



L 2 

c 

— ^ + C ||5u„| 
cn z 



L 2 



L 2 
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Since 



so 



and thus 



Ay 



< c 



L 2 



c- — u r 



L- 



T 't\ • ^(e|l«n|| ffi +§ + r 



n 



pi|<£ + (. K l H » + £ + ±) 



The term T 2 can be written as 



T 2 =(b(l- £ X ' + ) a (l - f A ' + ) ^-u n , i^j 



= T 2 1 +T 2 2 . 
Similar to the estimates for T 1; we have 



T 2 < 



and 



|T 2 2 | < C( 



1 

V'ej 



1>. 



ay 



Li, 



L 2 



<Cmax|a(£)| = 0(- 

lil>" 



a- — u 7 



L 2 



ci 1 



So 



Thus 



7 V r '( e KII„i+§ + ^ 



Re(&C A .P ej , (0C X u n ,u n 



< |Ti| + |T 2 



< + C ellu 
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Combining with (|3.8p . wc have 
Re ((A x - z) u n ,u n ) 



> { l-5)\ + C 5\\u n \t h -h -^-c(e»~ r ' + r " ' 



h ' ^ lrn|1 ^ 2 u c 2 \^ 1 n 2 ' n 

= i 5o -^ + (^-c £) |Kiy-c(^ + i 

> —5o, when n is large enough, 

by choosing e > and S £ (0, 1) such that e < ^-5 and <5 < ^SqIi. This finishes 

the proof of the lemma. I 

To study the essential spectrum of A x , we first look at the Zhislin Spectrum 
Z (A x ) ([IH])- A Zhislin sequence for A x and z £ C is a sequence {w n } £ H 1 , 
||w n || 2 = 1, supp u n C {£| |£| > n} and || („4 A — z) u„|| 2 — > as n — > oo. The set 
of all z such that a Zhislin sequence exists for A x and z is denoted by Z (A x ) . 
From the above definition and Lemma l3.5i we readily have 



Z(A X ) C jz £ C| Rcz > i(5 | • (3.11) 

Another related spectrum is the Weyl spectrum W (A x ) ((23)- A Weyl sequence 
for A x and z £ C is a sequence {u n } £ H 1 , \\u n \\ 2 = 1, u n — > weakly in L 2 
and || („4 A — z) it„|| 2 — > as n — > oo. The set W (-4 A ) is all z such that a Weyl 
sequence exists for A x and z. By (gHl Theorem 10.10]), W (A x ) C ct oss (A x ) 
and the boundary of a css (A x ) is contained in W (A x ) . So to prove Proposition 
ED it suffices to show that W (A x ) = Z (A x ). Since if this is true, then (|3.3p 
follows from (l3~TTi By ([29l Theorem 10.12]), the proof of W (A x ) = Z (,4 A ) 
can be reduced to prove the following lemma. 

Lemma 3.6 Given A > 0. Let \ £ Cq° (R) oe a cut-off function such that 
x\{\t\<B. } = 1 for some R > 0. Define \ d = x(£A0, d > 0. TTien /or 
eac/i d, („4 A — zj is compact for some z £ p (A x ) , and that there exists 
e (d) — > as d — > oo smc/i i/iai /or any u £ (R) , 

\\[A x ,Xd] u\\ 2 <e(d)(\\A x u\\ 2 + \\u\\ 2 ). (3.12) 

Proof. Since A x = J\f + IC X , where Af is positive and 

/C A = bC x P ey C x :L 2 ^L 2 (3.13) 

is bounded, so if z — —k with /c > sufficiently large, then z £ p (-4 A ). The 

compactness of Xd (A x + k) 1 follows from the local compactness of H 1 L 2 . 
To show (|3.12l) , we note that the graph norm of A x is equivalent to || • || H% . First , 
we write 



[lC X ,Xd} =b 



iXd 



P ey C A + bC A P e y 



nX 

L- ,Xd 
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We have 



C\xd 



1 - 



A 



X + V J tp, 



1 



-,Xd 



X 



x + v 



1 

^ey 



X 



x + v 



Xd,T> 



1 1 



X + Vtp. 



<•(/ 



—£ x ' + 
Xd 



< y 



Since 



is bounded, so 



C\ X d 



c 

< — 

l 2 ^l 2 Xd 



and therefore 



\[K\xd]u\\ 2 < 



C 



2 " Xd 



(3.14) 



Denote A/"i = 1 + ^ and M% is the Fourier multiplier operator with the symbol 



n 2 (k) 



tanh(M) (1 +ik)' 
Then TV = A/2A/1 and thus 

W, xd] = A/2 [M , X(i ] + [A/2 , xd] M . 

Since [A/"i,Xd] = ^x'(£/^) an d ||-^2|li2_»is i s bounded, we have 

C 

IIMEA/i.xJulla <-jNl a - 
To estimate [A/2, Xd], for w G Qf (R), we follow [HI p.127-128] to write 

[A/2, Xd] « = - (2^)~^ / n 2 (£ - y) (xd (0 - Xd (y)) » (y) dy 



(3.15) 



(2?r) 2 (f - y) n (£ - y) x d (p i€-y) + y)v (y) dydp 



^4 p t; dp, 



where A p is the integral operator with the kernel function 



K p (£, y) = - (2tt) Me - y) ^2 (£ - y) x'd (p (£ - y) + v) • 
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Note that a (£) = £fi2 (£) is the inverse Fourier transformation of in' 2 (k) and 
obviously n' 2 (k) 6 L 2 , so a (£) € L 2 . Thus 



// 



^ p (e,j/)| 2 dxd 2 / = 2^ y y |a| 2 (£ - y) 2 (p(C-y)+?/) dCd» 
= 2tt / / \a\ 2 (0\x'/(y) didy 



So 



and 



Thus 



= 27r|H|» a ||x , ji a = ^IHlLllX / ||^ 



ll[^2,Xd]|| L2 ^ L2 < 
|[A/- 2 , Xd ]A/"H| i2 < 5- ||Mu|| £a < 

C(2 6(2 

l|[Ar J x«j]«|| L2 <cfi + i > ) hi^. 



Combining above with (|3 . 14[> . we get the estimate (|3.12p . This finishes the proof 
of the lemma and thus Proposition [TJ I 

Recall that to find growing modes, we need to find A > such that A x has 
a nontrivial kernel. We use a continuity argument, by comparing the spectra of 
A x for A near and infinity. First, we study the case near infinity. 

Lemma 3.7 There exists A > 0, such that when X > A, A x has no eigenvalues 
in {z\ Rez < 0}. 

Proof. Suppose otherwise, then there exists a sequence {A rl } — * oo, and 
{k n } e C, {u n } 6 H 1 (R), such that Refc„ < and (A Xn — fc„) u n — 0. Since 
||.4 A — Af\\ = ||^ A || < M for some constant M independent of A and TV is a 
self-adjoint positive operator, all discrete eigenvalues of A x lie in 

D M = {z\ Rez> -M and \Imz\ < M} . 

Therefore, k n — * koa G Dm with Refcoo < 0. Denote 

e(0=nua{|o(0|,|6(0|. |5(0l}, (3-16) 

where (£)>&(£)> c(£) are defined in and (fXTUj) . Then e (£) -> as 

|£| — » oo. Define the e (£)-weighted L 2 space L 2 with the norm 

\\u\\ Li = (Je(0\u\ 2 deV. (3.17) 

We normalize it„ by setting ||u n || i2 = 1. We claim that 

||w n || i < C, for a constant C independent of n. (3.18) 

H 5 
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Assuming (I3.18p . we have u n — > weakly in H 2 . Moreover, Uoo 7^ 0. To show 
this, we choose R > large enough such that ma,x^ >R e (£) < Then 

1 



Since 



2 1 
|£|>fl zo 
strongly in £ 2 ({|£| < i?}), we have 



< 



t\<R 



• (£) |iioo| 2 = lim 



ISI<« 



e(£)K| 2 ^> \ 



and thus Woo 7^ 0. By Lemma I3~2| A Xn — > A/" strongly in L 2 , therefore .A An ii n — * 
Muoo weakly. Thus Muoo = kooUoo. Since Refcoo < 0, this is a contradiction 
to that A/" > 0. It remains to show (|3.18p . The proof is quite similar to that of 
Lemma I3~5l so we only sketch it. From (_4 An — k n ) u n = 0, we have 



(Afu n ,u n ) + Re(bC x -P ey (0C Xn u n ,u n ) = Rek n \\u n \\ 2 2 <0 



By Lemma T3.1[ 



(Afu n , Un) >(1-S)- \\u n \\ 2 L 2 + C Q S ||it n ||^i 



Following the proof of Lemma 13.51 we write 
bC x "P ey (£)C x "u n ,u 

9 ( Hey 



1 



Hey ( 



1_£A„,+ 



4>e 



bip e y a, 1 - £ 

_ Tl + T? + T% +T$, 
The first term is estimated as 
1 



< u 



l-£ x '» 

l-£ Xn '- 
1 — £ Xn ' 
(l-£ x 



t'J/ 



;j [l\\u n \\^+C\\u n \\ 2 L l 



b- I Ti rl 1 2 dx 

Wey 



- \\u n \\\* + C IIm^II^ 



< 9 — Il«n|li2 + C||u„|| i2 — ||u„|| L2 +C||u n | 



C2 

< ^ \\UnWl2 +C||U„|| L2 \\u n \\ L 2 



< — \\ u n\\ L 2 + £ \\Un\\ L i + C e ||u„|| i2 



C\\u n \ 



L 2 



(3.19) 
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where in the second inequality, we use the fact that 



The second term is controlled by 
I?? | <C( t(l-£ x ^+' 



1 1 




1 1 


-0 ey 3 C 3 




tpey C 



<Ce(0- 



+ \\ u n\\ L 2 ) \\Un\\ L 2 



2 

LI ' 



< e ||n„|| i2 + C s \\u n \\ L2 



< C (e\\u n \\ H i + \\u n \\ L2 ^j \\u n \\ L2 < Ce\\u n \\ 2 H i + C £ \\u n \\ 
where in the second inequality we use Lemma 13.41 (iv). The third term is 

1^2 | < C\\u n \\ L 2 \\U n \\ L 

By the same estimate as that of Tj 2 , we have 

|T 2 2 | < Ce\\u n \\ 2 Hh +C e \\u n \\ 2 Ll . 
Plugging all of the above estimates into ()3.19|) . we have 

> 



<'-*s-£ 



\\u n \\ 2 L 2 + (C S - Ce) \\u n f H i - C e ||m„||^2 



> 7^0 IKIIi2 + ^C S \\u n f H i - C e ||w„||ia 



by choosing 6, e such that 

c 1 rjr 1 C S 

Then ([3~T8| follows. ■ 

4 Asymptotic perturbations near zero 

In this Section, we study the eigenvalues of operator A x when A is very small. 
By Lemma [321 when A — > 0+, A x — > .4° strongly, where 

The related operator in the physical space is .4° : i? 1 (<S e ) — > L 2 (S e ) defined by 

which is the strong limit of _4g when A — > 0+. We have the following properties 
of .4. . We use A (/J.) to denote the dependence on the solitary wave parameter 
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Lemma 4.1 (i) The operator A : H 1 (R) — * L 2 (R) is self-adjoint and 

(ii) ift ex (£) € ker„4° and A has at least one negative eigenvalue that is 
simple. 

(iii) Under the hypothesis (HI) of no secondary bifurcation, k.erA°(p) — 
{ip ex (£)} when /i is not a turning point and ker*4° (jjl) = ^ip ex (£) ; } 
when n is a turning point. For any /i > ^, -0 e:c (£) is the only odd kernel of 
A°(fx). 

(iv) When fx — — is small enough, A (fi) has exactly one negative eigenvalue 
and ker*4° (/i) = {ip ex (£)}. Under hypothesis (HI), the same is true for A a (/i) 
with /i G where is the first turning point. 

(v) When fi — > oo, </ie number of negative eigenvalues of A (/i) increases 
without bound. 

Proof, (i) The essential spectrum bound follows from the observations that 
Cess (A/") = [^,+00) and bP ey /ipl y -> - Jr when |£| -> 00. 

Proof of (ii): To show tp ex (£) € ker„4°, it is equivalent to show that 

ipex (x) = ipex {x,rj e (x)) € ker^4°. 

On S e , we have 

'ipex (x) + r]ex^ey (x) = 0, P ex (x) + ?7 ex P ey (x) = (4.1) 

and 

Pex \X) — {^Pex^Pexx ~t~ ^ey^Peyx) — ^Pey {j]ex^eyy ^eyx) (^*2) 

So 

^-Veo: (x) = -^r-Vex = ^ = (<Pex) = -Me (ipex (x)) , 

^ej/ Vey Vey «X 

and thus ylgV'ex (x) = 0. Now we show that .4° has a negative eigenvalue. We 
note that the Fourier multiplier operator J\[— h has the same symbol as in the 
Intermediate Long Wave equation (IIW), for which it was shown in [2] that for 
K > large, the operator (J\f+K)~ is positivity preserving. Thus, by the 
spectrum theory for positivity preserving operators ([2]), the lowest eigenvalue 
of A is simple with the corresponding eigenfunction of one sign. Since ip ex (x) 
is odd, ip ex (£,) has a zero at £ = 0. So is not the lowest eigenvalue of .4° and 
.A has at least one simple negative eigenvalue. 

To prove (iii)-(iv), first we show that the operator A (fi) is exactly the oper- 
ator A (A) introduced by Plotnikov ([581 P- 349]) in the study of the bifurcation 
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of solitary waves. In [S5J, h is set to 1 and the parameter A = is the 

inverse square of the Froude number, then the operator A (A) is defined by 

A (A) = Af—a, a (?) = A exp (3r) cos 6 + 6' (?) 

In the above, exp (r + i#) = where is defined by (|2. 18[> . as can be seen 
from ( 58, (4.2), p. 348]) with u = w. To show that A (fi) = A(X (/x)), it 
suffices to prove that 

A (A) tjj ex (£) = 0. (4.3) 

Since this implies that 

= (A (A) - .4° (/i)) ^ (0 = {-a - bP ey /^ 2 ey ) $ ex (0 

and thus bP ey /ip 2 y = —a. We prove (|4.3[) below. In 58J, solitary waves are 
shown to be critical points of the functional 



J{\w) = \l {wNw-Xw 2 {l+Nw)}d£. (4.4) 
2 /r 



Let the self-adjoint operator Aq (A) to be the second derivative of J (A, w) at a 
solitary wave solution. In [55J P- 349], the operator A (A) is defined via 

A (A) = M*A (A) M. 

Here, the operator M : L 2 — > L 2 is defined by 

Mf = / (1 + Cuj') + w'C/ = Re {Wft/} , (4.5) 

where C is defined in Section 2 such that IZf — f — iCf is the boundary value 
on {<; = 0} of an analytic function on Dq. Our definition (|4.5[) above adapts the 
notations in [T3J p. 228] , which studies the bifurcation of Stokes waves by using 
a similar variational setting as [S5J. Taking d/d£ of the equation V W J (A, w) = 
for a solitary wave solution w, we have Aq (A) u/ = 0. Since 

M V = Re <^ ^ = Re ^ 

( J \l + Cw' + iv 



Re- 



(u/ - iCw') (1+Cw' - w') 

I 1 

,2 = ~ V * = --Wex (?) , 



we have Mip ex (?) = — cw' and thus 

^4 (A) Vex (0 = -cM*i4o(A)«/ = 0. 
This finishes the proof that A (p,) = A (A). 
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Proof of (iii): By applying the analytic bifurcation theory in [TJ], [T5] to the 
variational setting (|4.4[) for the solitary waves, one can relate the secondary bi- 
furcation of solitary waves with the null space of A (equivalently V^J) .Under 
the hypothesis (HI), there is no secondary bifurcation and therefore the kernel 
of .4° is either due to the trivial translation symmetry (^ ex ) or due to the loss 
of monotonicity of A (p) at a turning point which generates an additional kernel 
9 M Ve- In the Appendix, we prove that at a turning point py, A d^ipe = 0. By 
|23j there is no asymmetric bifurcation for solitary waves with F > 1, so i\) ex (£) 
is the only odd kernel of ^4° (p). 

Proof of (iv): Let A (p) = exp (— 3p 2 ), then p rj ^ is equivalent to A » 1 
and thus p is a small parameter. Consider an eigenvalue v of A (/i(A)), let 
v = p 2 (3 — a (p)). By using the KDV scaling, it was shown in [58j p. 353] that 
when p — > 0, the limit a (0) is an eigenvalue of the operator 

which has three eigenvalues §, 3 and 2Z. Therefore, when p is small, .4° has three 
eigenvalues — ^-p 2 + o (p 2 ) , o (p 2 ) and \p 2 + o (p 2 ) . Since is an eigenvalue of 
A , the middle one must be zero and the rest two eigenvalues are one positive 
and one negative. Under the hypothesis (HI), when p < pi, that is, before the 
first turning point, we have ker^4° (p) = {ip ex (£)}. Then for all p € (^, p\) , the 
operator ^4° (p) always has only one negative eigenvalue. Suppose otherwise, 
then when p increases from - to p\, the eigenvalues of A (p) must go across zero 
at some p = p* € (-, pi) .This implies that dimkerA (p*) > 2, a contradiction 
to (HI). 

Property (v) is Theorem 4.3 in [58]. I 

We note that by Lemma [4~7T1 (iv), there is no secondary bifurcation for small 
solitary waves. Although this fact was not stated explicitly in [55], it comes as 
a corollary of results there. 

Next, we study the eigenvalues of A x for small A. Since the convergence of 
A x — > .4° is rather weak, we cannot use the regular perturbation theory. We 
use the asymptotic perturbation theory developed by Vock and Hunziker ( [67]). 
see also [25], [30] . First, we establish some preliminary lemmas. 

Lemma 4.2 Given F S Qj° (R). Consider any sequence X n — > 0+ and {u n } € 
H 1 (R) satisfying 

||^ A "u„|| 2 + ||u„|| 2 < Mi < oo (4.6) 
for some constant Mi . Then if w — linin^oo u n — 0, we have 

lim ||Fu n || 2 = (4.7) 

n — ^oo 

and 

lim ||U A ",Flu„|| =0. (4.8) 
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Proof. Since (|4.6[) implies that < C, (|4.7[) follows from the local 

compactness of H 1 L 2 . For the proof of (|4.8p . we use the same notations as 
in the proof of Lemma ETUI We write A Xn = TV + /C A " . Then 

[A/ - , F] - [A/2M , F] = A/2 [M , F] + [A/2 , F] m , 

where A/1 = 1 + ^ and A/2 has the symbol ri2 (fc) defined by (|3.15p . We 
have || [A/1, F] u n || 2 = ||F'u n || 2 — + 0, again by the local compactness. Since 
Ju"2 (k) — > when |fc| — > 00, by [HI Theorem C] the commutator [A/2,F] : 
L 2 — > L 2 is compact. This can also be seen from the proof of Lemma 13.61 
since [A/a > F] = f Q A p dp where A p is an integral operator with a L 2 (R x R) 
kernel. Since ||A/iu n || 2 < ll u n||#i < C and u„ — > weakly in L , we have 
w„ = A/l^n -> weakly in F 2 . Therefore, || [A/2, F] A/1m„|| 2 — | [A/2, F] u„|| 2 — > 0. 
Thus, ||[A/\F]u n || 2 -> 0. Since 



[/C A ",F] 



( 1 -*" + )&( 1 -* x " + ) 



F£ 



A„,+ 



fpey 

1 _ ( ?A„,+ 



"0ey '" 



1 



6 1-5 



A„,+ 



F, 



F,£ 



A„,+ 



= Pn + Qn- 

For any e > 0, by Lemma 13.41 (iii), when n is large we have 



F£ 



l2 <e|K||^ <e||«n|| H i <eAfi. 



So 



\1n\\ 2 



< C . 



F£ 



\Fu. 



nil 2 



< eCAF + C ||FuJ| 2 < 2eCM u 



when n is large. By the same proof as that of (|3.5p . for any A > 0, we have the 
estimate 



?a,+ 



Denote 



iPey V 

then II^H^i < C||m„|| h i < CM\. Since w — lim. 
'„ = as 

of 1 1 C?n. 1 1 2 ) we nave 



< C (independent of A) . 



1 



U n , 



n — >oo 



0, we have w 



Iim n _ >ac r n = as in the proof of Lemma l3.4l (iii). Then similar to the estimate 



IIP* 



F£ 



|Fr„ 



< 2eCM u 



when n is large. Therefore, || [/C A ",F] u n \\ 2 < AeCMi when n is large enough. 
Since e is arbitrary, we have || [/C A ",F] u n || 2 — > 0, when n — > 00. This finishes 
the proof of gH). ■ 



Lemma 4.3 Let z G C with Rez < *8o, then for some n > and all u G 
Co° (\ x \ — n )> we have 

\\(A X -z)u\\ 2 >±S \\u\\ 2 , (4.9) 

when A is sufficiently small. Here 5q > is defined by Jff.^l ). 
Proof. The estimate f|4.9[) follows from 

Re((A x - z)u,u) > ~5 \\u\\l. (4.10) 

The proof of (14. 10[) is almost the same as that of Lemma l3.51 except that Lemma 
13.41 (iii) is used in the estimates. So we skip it. I 

With above two lemmas, we can use the asymptotic perturbation theory 
( 29 , 30 ) to get the following result on eigenvalue perturbations of A x . 

Proposition 2 Each discrete eigenvalue ko of A with ko < iSo is stable with 
respect to the family A x in the following sense: there exists Ai, 6 > 0, such that 
for < A < Ai, we have 
(*) 

B{k ;S) = {z\ < |*-*o| < 5} CP(A X ) , 

where 

P (-4 A ) = |z| R x (z) = (A x — z^j exists and is uniformly bounded for A G (0, Ai) 
(ii) Denote 

P x = f R x (z) dz and P = f R° (z) dz 

J{\z~k \=S} J{\ z -k \=S} 

to be the perturbed and unperturbed spectral projection. Then dim P\ — dim Pq 
and liniA^o \\P\ - Poll =0. 

It follows from above that for A small, the operators A x have discrete eigen- 
values inside B (ko; S) with the total algebraic multiplicity equal to that of ko. 

5 Moving kernel and proof of main results 



To study growing modes, we need to understand how the zero eigenvalue of 
A is perturbed, in particular its moving direction. In this Section, we derive 
a moving kernel formula and use it to prove the main results. We assume 
hypothesis (HI) and that fj, is not a turning point. Then by Lemma [4.11 (iii), 
ker^4° (fi) — {ip ex (£)}. Let Xi,5 > be as given in Proposition [5] for ko = 0. 
Since dim P\ = dim Po = 1, when A < Ai there is only one real eigenvalue of A x 
inside B(0;S), which we denote by k\ G R. The following lemma determines 
the sign of k\ when A is sufficiently small. 
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Lemma 5.1 Assume hypothesis (HI) and that /i is not a turning point. For 
A > small enough, let k\ G R to be the eigenvalue of A x near zero. Then 

where E (/i) is the total energy defined in il.8\) . 

The following a priori estimate is used in the proof. 

Lemma 5.2 For A > small enough, consider u G H 1 (R) satisfying the equa- 
tion (A x — z) u — v, where z G C with Hez < ^So and v G L 2 . Then we 
have 

Ny <c(\\u\\ L2c + \\v\\ L2 ), (5.2) 

for some constant C independent of X. Here, the norm \\-\\ L2 is defined in \3. 1 7\ ) 
with the weight e (£) defined by \3.10jl . 

Proof. The proof is almost the same as that of the estimate p. 181) in the 
proof of Lemma [377] So we only sketch it. We have 



(Afu,u) +Re (bC x P ey (£)C A u,u) - Rez||u||? = Re (u, v) . 

By the same estimates as in proving (|3.18|) . except that Lemma l3~4l (iii) is used, 
we have 



||u||£ a +(Co6-Ce)\\u\\ 2 H i-C £ \\u\\ 2 Ll 



c-"s-^-K 

<e\\u\\ 2 L2 + - \\v\\ 2 L2 . 

Then the estimate (|5.2I) follows by choosing e > and S G (0, 1) properly. I 
Assuming Lemma I5.R we prove Theorem [TJ 

Proof of Theorem [l] We fixed /i G (fit, fix). Under the assumption 
(HI), it follows from Lemma H. II that A (fi) has only one negative eigenvalue 
Icq < and ker^°(^) = {tpex (&}■ By Proposition [5] and Lemma fSTTT there 
exists Ai,<5 > small enough, such that for < A < Ai, A x has one negative 
eigenvalue in B (k^ '; 5) with multiplicity 1 and one positive eigenvalue k\ in 
B (0; 5) because = E' (fx) jd (/x) < for /i G (jtti, fix). Consider the region 

Q, = {z\ > Rcz > -2M and |Imz| < 2M}, 

where M is the uniform bound of ||.4 A — 7V|| . We claim that: for A small enough, 
A x has exactly 2 eigenvalues (counting multiplicity) in 

D.S = {z\ 28 > Rez > -2M and \Imz\ < 2M} . 

That is, all eigenvalues of A x with real parts no greater than 28 lie in B [k^ ; 8) U 
B (0; 8). Suppose otherwise, there exists a sequence A„ — > 0+ and 

{«„} G H 1 (R) , z„ G £1/ (B (ko;S) U B (0; 6)) 
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such that (A Xn — z n ) u n = 0. We normalize u n by setting ||ii n ||z,2 = 1. Then 
by Lemma l5.21 we have ||u n || 1 < C. By the same argument as in the proof of 
Lemma |3. 71 u n — > Uoo ^ weakly in H s. Let 

z n ^ Zoo en/ (B (k ;S) UB(0;S)) , 

then A°u OQ — ZaaUao which is a contradiction. This proves the claim. Thus for 
A small enough, A x has exactly one eigenvalue in fl. 

Suppose the conclusion of Theorem [1] does not hold, then A x ([i) has no 
kernel for any A > 0. Define nn (A) to be the number of eigenvalues (counting 
multiplicity) of A x in ft. By (|3.3p . the region tt is away from the essential 
spectrum of A x , so nn (A) is a finite integer. For A small enough, we have 
proved that nn (A) = 1. By Lemma \3. 71 nn (A) = for A > A. Define the two 
sets 

Sodd = {A > 0| n n (A) is odd} and Seven = {A > 0| iifi (A) is even} . 

Then both sets are non-empty. Below, we show that both S dd and Seven are 
open. Let Ao £ Sodd and denote fci,--- ,ki (I < nn (Ao)) to be all distinct 
eigenvalues of A x ° in fl. Denote ihi, • • • , ih m to be all eigenvalues of A x " on 
the imaginary axis. Then \hj\ < M, 1 < j < m. Choose S > sufficiently small 
such that the disks B (ki, S) (1 < i < I) and B (ihj; S) (1 < j < m) are disjoint, 
B(ki\5) C £1 and B(ihj\5) does not contain 0. Note that A x depends on A 
analytically in (0, +00). By the analytic perturbation theory ([29]), if |A — Aq| is 
sufficiently small, any eigenvalue of A x in lies in one of the disks B (fc^; 5) or 
B (ihj; 5). So nn (A) is nn (Ao) plus the number of eigenvalues in U^ 1 S (ihj; 5) 
with the negative real part. The later number must be even, since the complex 
eigenvalues of A x appears in conjugate pairs. Thus, nn (A) is odd for |A — Ao| 
small enough. This shows that S^dd is open. For the same reason, S^von is 
open. Thus, (0, +00) is the union of two non-empty, disjoint open sets Sodd and 
Seven- A contradiction. 

So there exists A > and / 11 £ H 1 (R) such that A x u = 0. Let / = 
B^u e H 1 (S e ), then A x f = 0. Define n{x) = C x f, P{x) = -P ey n(x) and 
ip (x, y) to be the solution of the Dirichlet problem 

At/> = in V e , 7/>| 5e = / , ^(.t, -h) = 0. 

Then (n (x) , -0 [x, yj) satisfies the system (|2.3j) - (|2.7[) . thus e xt [n (x) , ip (x, y)] is 
a growing mode solution to the linearized problem (|2.2| . Below we prove the 
regularity of [77 (x) , ip (x, y)] . Since the operator C x is regularity preserving, from 
A x f = we have 

Mx) = ~C x P ey C x f eH 1 (S e ). 

By the elliptic regularity of Neumann problems ([!]), we have ip (x, y) S H 5 / 2 (V e ) 
So by the trace theorem, / = V|s c e H 2 («^e)< Repeating this process, we get 
ipn{x) S H 2 (S e ) and ip (x,y) € H' 7 1 2 (T> e ) . Since the irrotational solitary wave 
profile and the boundary S e are analytic ( [36] ) , we can repeat the above pro- 
cess to show ip (x,y) £ H k (V e ) for any k > 0. Therefore n (x) — C x (ip\s e ) & 
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H k (S e ), for any k > 0. By Sobolev embedding, [rj (x) , ip (x, y)] € C°°. This 
finishes the proof of Theorem [1] I 

Proof of Theorem [2j Let /xi < /Z2, • • • < /Lt„ < • • • be all the turning 
points. Then \i n — ► +00. Under the assumption (HI), ker.4 (fa = {^ez (£)}, 
for /i 6 (jUj, /ij+i) , i > 1. Denote by n~ (/z) the number of negative eigenvalues 
of .4° (/i). Then n~ (fa is a constant in (/Xj,/Xi+i), by the same argument in 
the proof of Lemma 14.11 (iv) . Denote fa < fa, ■ • • < fa < ■ " to be all the 
critical points of E (/i). Each fa- lies in some interval (fj,i, /Uj+i). Then the sign 
of dE/dc = E' (fa jc' (fa changes at fa in (/ij, /Ui+i). So we can find an interval 
It C (fa, (J-i+x) such that the number 

h- (fa = n- (m) + (1 + si<?n (M) M (M))) /2 

is odd for Note that n~ (//) is the number of eigenvalues of A x (fa in the 

left half plane, for A sufficiently small. So by the same proof as that of Theorem 
[TJ we get a purely growing mode for solitary waves with [i S Ik- Since fa — > 00, 
the intervals I n goes to infinity. I 

We make two remarks about the unstable solitary waves proved above. 

Remark 1 In terms of the parameter u> = 1 — F 2 q^, it was found (164-1) from 
numerical computations that the energy maximum is achieved atujtx 0.88, which 
corresponds to the amplitude-to- depth ratio a = rj e (0) /h — 0.7824 ('521). The 
highest wave has the parameters u = 1, a = 0.8332 and the maximal travel speed 
is achieved at to = 0.917, a = 0.790 ( \bl$ , JT^jj)- So the unstable waves proved 
in Theorem [7] and fj| are of large amplitude, and their height is comparable to 
the water depth. Therefore, this type of instability cannot be captured in the 
approximate models based on small amplitude assumptions. Indeed, although 
some approximate models share certain features of the full water wave model, no 
unstable solitary waves have been found. For example, the Green-Naghdi model 
is proposed to model water waves of larger amplitude and it also has an indefinite 
energy functional, but the numerical computation ( \31\ P. 529] ) indicates that all 
the G-N solitary waves are spectrally stable. For Camassa-Holm and Degasperis- 
Procesi equations, there exist cornered solitary waves (peakons). However, these 
peakons are shown to be nonlinearly stable (JEEl, 14°V )- So the instability of large 
solitary waves seems to be a particular feature of the full water wave model. 



Remark 2 The linear instability suggests that the solitary wave cannot preserve 
its shape for all the time. The long time evolution around an unstable wave was 
studied numerically in ]65jj . It was found that small perturbations with the same 
amplitude but opposite signs can lead to totally different long time behaviors. For 
one sign, the perturbed wave breaks quickly and for the other sign, the perturbed 
wave never breaks and it finally approaches a slightly lower stable solitary wave 
with almost the same energy. Note that in the breaking case, the initial perturbed 
profile has a rather negative slope ([65]). The wave breaking for shallow water 
waves models such as Camassa-Holm (1 21]/ ) and Whitham equations ([60], \51$ ) 
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is due to the initial large negative slope. It would be interesting to clarify whether 
or not the wave breaking found in J65]j has the same mechanism. 

The wave breaking due to the instability of large solitary waves had been 
used to explain the breaking waves approaching beaches fJ25f. J56"jj . 157)1 ). When 
a wave approaches the beach, the amplitude-to- depth ratio can increase to be near 
the critical ratio (m 0.7824,) for instability and consequently the wave breaking 
can occur. 

It remains to prove the moving kernel formula (|5.1[) . 

Proof of Lemma 15.11 As described at the beginning of this Section, for 
A > small enough, there exists u\ € H 1 (R), such that (A x — kx) u\ — with 
k\ € R and liniA^o+ kx = 0. We normalize u\ by ||ua|| ra = !• Then by Lemma 
15.21 we have ||ua|| i < C and as in the proof of Lcmma l3.7[ u\ — > uq ^ weakly 

H 2 

in H?, Since A°uo = and ker*4° (/x) = {tp ex (£)}, we have uo = CQip ex (£) for 
some Co ^ 0. Moreover, we have \\u\ — uq\\ i = 0. To show this, first we note 
that ||ma — uo|| L 2 — > 0, since 

||tfrA — i*o||i» < / e (£) \u\ - u \ 2 d£ + max e (£) ||u A - u \\ 2 L2 , 
./|f|<R ICI>-R 

and the second term is arbitrarily small for large R while the first term tends 
to zero by the local compactness. Since 

(.4 A - k x ) (ux - u ) = k x u + (A - A x ) u , 

by Lemma 15.21 we have 

IK - uo|| H i < C (||u A - uollij + N HuAllia + || (-4° - A x ) u \\ 2 L2 ) -» 0, 

when A — > 0+. We can set Co = 1 by renormalizing the sequence. 

Next, we show that limA^o+ % = 0. From [A x — kx) u\ = 0, we have 

u x A x -A° k x 
A — H u\ = (5.3) 

Taking the inner product of above with if> ex (£), we have 
k x /A x -A° 



x (ux,1pex{0)= ( ^ U>\,4>ex(£)J '■= TTl (A) . 

We compute the integral m (A) in the physical space, by the change of variable 
£ — > x. Denote by ( , ) the inner product in L 2 (S e ). Noting that dx = b (£) d!;, 
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we have 



' A x - A 
m W = ( — ^ — ~ u a (a;) , V'ex (a:] 



1 Pe 9 , \ Pey 1 1 

X + V^j y I \lpey X + V 1p ey 

1 Pey c \ J- 1 , 



X + V 1p ey 1p ey 

= mi + mi + "la, 

where we use 

(<•»,.. W C»-£jj) (5.4) 



nj 



1 p p 1 1 IP 1 

1 r ey _ r ey ± -I x ^ey f A,+ 1 



A + I? l/'ey V'ey A + V lp ey X + T> 1p ey ~ lf) ey 

We compute each term separately. For the first term, when A — ► 0+, 

mi (A) = (xT^W v Ux ^ eyr]ex ) = (&" A ' ~x~i5^ ex 

Wey Wey I \Wey Wey / 

where we use Lemma 12.11 (b) in the above and the resultant integral is zero 
because P ey ,ip ey ,r] e are even and ip ex is odd in x. The second term is 

1712 (A) = aT^ Ua '^) = ( XT^^'^i ( ^ } 

-> - (Wex, -J— Ue ) = 0, 
\ Vew / 

where the relations (|4.ip and (|4.2p are used in the above computation. For the 
last term, 

1713 (A) = - ( xvd^ 8 ^^^) = - ( Pe ^ + jr v ux > x^v 71 -) 



P ey £ x '+^-u x , (£ x - - 1) -L Vl 



0, 
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because E x ^ -r—u\ — > weakly in L 2 , when A — ► 0+. So m (A) — » 0, and thus 

v kx r TO ( A ) n 
hm — = lim ■ -— = 0. 

A^0+ A A^0+ (u X ,1pex (£)) 

Now we write u\ = c A -0ei + A«a, with c A = {ux,ipex) / (ipex,i>ex)- Then 
(v\,ip ex ) = and ca — » 1 as A — > 0+. We claim that: ||f a||^2 < C (independent 
of A). Suppose otherwise, there exists a sequence A„ — » 0+ such that ||ua„ 11^2 > 
n. Denote v\ n — v\ n j \\v\ n \ L i- Then \v\ n || i2 = 1 and V\ n satisfies 

1 fk\ A x ™ -A \ 

A x -dx n = ^u Xn - c Xn " ^ ex (0 . (5.5) 



Denote 



where 



I«aJIl2 



A x ™ - A 

9n (0 = r (0 = & (0 »a„ (a; (£)) . 

A n 



A x - A 

W\ (x) = — ^ -^ex (x) (5.6) 

J -Pea Pey 1 j_ _, C A.+ jefl 1 1 \ , 

./. \ i t-i ./. I Vey'lex 



A + 2? Veu i'ey ^ + T> 1p ey 1p ey \ + T> %p ey 



1 ^ / \ , ^e-y f r-A 4- 1 



ex 



A + 2? cte -0e a A + 2? Vet, A + 2? 

(1 _ JffL + £k (i _ Jfe_ _ £ x ,+ ^v (i _ £ A,+) .^L 



By Lemma [2TT1 ||i0a||x,3(5 ) < C (independent of A), and moreover 

w\ (x) — > — — H L -^- strongly in L 2 {S e ), when A — > + . (5.7) 

So Hffnllia < C and thus by applying the estimate (|5.2|) to (|5.5p . we have 
||wa„|| i < C. Therefore, as before, v\ n ^ io ^ weakly in H?. Since 
"T*S ||- 1 ii > 0, we have yl «o = 0. So uo = cii/i e2; (£) for some ci ^ 0. But 

(ibi„, V'ex (£)) = implies that (vo,ip e x (£,)) = 0, a contradiction. This proves 
that ||i>a||£2 < C. The equation satisfied by v\ is 

A X V X = ^U X - CA ^ ^ Ipex (0 = T^-WA - C A fo (0 WA (£)) • 
A n A A n 

Applying Lemma [5.21 to the above equation, we have \\v\\\ i < C and thus 



ua — ► weakly in 22 2 . By (|5.7p . Uq satisfies 
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It is shown in the Appendix that 

A%j e (x) = -(^ + ^)(x), (5.8) 
\Vey Vtu J 



>Jey 

and equivalently, 

\Vey Vey J 

where 

i) e (x, y) = ip e (x, y) - cy. 

So A (y - d c ij e (£)) = 0. Since (v , ip ex (£)) = lim A ^ + (v\, Vex (£)) = 0, we 
have 

V = d c tp e (£) + d tp ex (£) , d Q = - (d c ^ e , 1p ex ) I \i>em\i,i ■ 

By the same argument as in the proof of IIua — Wo II i — » 0, we have ||u A — foil i 
0. We rewrite 

"A = Cxlpex + AWA = C\1p ex + \v\, 

where c\ = c A + Xd , v\ — v\ d ip ex . Then c\ — > 1, v\ — > <9 c Vi e (£) when 
A — > 0+. 

Now we compute limA_»o+ From (|5.3p . we have 

u x A x -A° fc x \ fc A 

A^" A \~A~ x / = A2 UA ' 

Taking the inner product of above with i\) ex (£) , we have 

{U\,1pex {£)) = CA I T2 Vex, Vex I + I T «A, % 

= C A /l + I 2 . 

Again, we do the computations in the physical space. For the first term, we use 
CT) to get 

/ A x -A° e \ / w x (x) 

h = ( ^ ^ W ' W / = \ ^ > ^ W 

/ V U e L , / Pey 2? Ve 



Vex (x)) + ( - — - — ,^ ei (a;) 



\ (A + P) A 7/; ey ' re " v 7 \> e „ (A + X>)AVe 
= I{ + J? + 
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We have 

/i1 = -((x-aTp)^'^) 
11 = {~A V (a ~ \Tv) it v ,iJevVex ) = { Pey \ + v^ ex ) 

= (xr^ity^ ey i {Ue) ) = {Ue) ) = (*> - l) it v 

( ??e ' Ay)' 

and 

= (1 " ^ A~y' {£X ^ ~ l)J t v )^- ( PeV ity' A~) ' 

So 

lim h (A) = - ( T) e , + ) • 

A^0+ \ Ay / 

To compute I 2l we write 

r / 1 p ey - , \ , 1 1 _ . \ 

J 2 = \ \ i ^ -T V ^^eyVex ) + ( J 1 , , ^ -p- V X , Ay Vex > 

~ \ \~~^rr LS X,+ -^—V X ,1peyVex ) 
XA + I?^ Vey / 

= jl+jf +/f. 

We have 

J » = t^*-) = (t^< (£ x - - l) ^\ - - (^fl c fc, J5-N 

\ ^ A - X» / \ Ay Ay/ \ Ay Ay I 

^•<- 1+ ^£>--<**'£> 
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and 



r 3 - - 

1 2 — 



1 



< y 



Pey£ 



1 



1>. 



-v\, 



■i + e 



'■y 



A,-\ Ve_ 
Ipey 



So 



Thus 



lim I 2 (A) = - ( d c Tp e , -P- + 



p. 



,. k x cxh + h 

lim -rr = lim 



A^0+ A 2 A^0+ (ux^ex) 



n 



Pey 



ex || L2 



It is shown in the Appendix that 



u, 



P 



where the momentum P is defined by 



2 

ej/ 



(5.9) 



ex. It is shown in [10] (see also [47]) that for a solitary wave 

dP 



with <p e — (f. 
solution 

dE _ 
dc dc ' 

where we note that the travel direction considered in this paper is opposite 
to the one in the above references. A combination of above results yields the 
formula (j5Jj) . I 

As a corollary of the above proof, we show Theorem [3J 

Proof of Theorem [3[ The proof is very similar to that of Lemma 15.11 
so we only sketch it. The main difference is that the computations depend on 
the parameter \i n . We use rj e , n ,t/j e y, n etc. to denote the dependence on [i n , 
and rj e ,ip ey etc. for quantities depending on fig. Denote u n (£) = B n {ip n \s e „)• 
Then A„ n u n = and we normalize u n by ||u„|| i 2 =1- First, we show that 

£^™' ± — * strongly in L? (R). Indeed, for any v € L 2 (R), since the norms 



I 'Ilia 



and 



are equivalent, we have 



< c 



c 



A 



c 



d\\M a , n v 



-d\\M a , n v\\ 2 L 2 



a\<S 



Ai + A 2 , 
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where 6 > is arbitrary and ^M a n :a € R 1 j is the spectral measure of the 

self-adjoint operator R n — —iV n on L\ n ^, iy n - For any e > 0, since Af{o} = we 
can choose S small enough such that 



d\\M a vf Ll <§. 



Since the measure d||M Qjn ?;||^ 2 converges to d|]M Q w||^ 2 , we have A 2 < 

h n^ey,n ^ &y 

e when n is big enough. The term A\ tends to zero when n — > oo. Because e 
can be arbitrarily small, we have 



lim 



2 

= 0. 

L 2 



So A* n — > -4° strongly and similar to the proof of Lemma 15.11 we have u n — > 
tpex (0 in -ff 5 (R) by a renormalization. We write u„ = c n ip ex ^ n (£) + A n i; n ,with 
c n = (u n , ipex,n) / {4>ex,n, ^ex,n) ■ Then c„ -> 1 and {v n >i>ex, n ) = 0. As before, 
we have v n —> v in H ? (R), where 

V = d c t[) e (£) + doV'ex (0 , ^0 = ~ (c?c V^e , ^ex) / (4>ex,1pex) ■ 

We rewrite 

^ra C-nlpex.n is) ~\~ ^n^n — Cn^ex.n ~t~ A n V n , 

where c n = c„ + X n d and W n = v n - doip e x,n- Then c„ — > 1, W„ — > d c 4> e (£) when 
n — ► oo. Similarly, we have 

A\n _ /10 \ / /(A„ _ 40 

■^Vi <^ n 1 1 1 1 f n 

ex ,n -i y ex .n I i I r 
A„ / \ A„ 



= C„ ( -1pex,n,fpex,n ) + ( — ^ ~V n , 1pex,n ) = <kih + h- 



As before, the calculations of I\, 1% are first done in the physical space S e ^ n with 
the inner product (., .) . We use the same notations as in the proof of Lemma 
15.11 By the same computations, we have 

\ tyey,n / n \ ^ ' Wey.n / 

Vey J \ Vey 

The other terms are handled similarly. So finally we have 

= lim C n h +h=-{ Ve,2— + -pr'He ) - { O c 1pe, h -pf^ 

n ^°° \ Vey Vey I \ Wey Vey 

IdE 

"7 Mo ■ 

c dc 

So E' (mo) = 0. This finishes the proof of Theorem H I 
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6 Appendix 



In this Appendix, we prove (|5.8[) and (|5.9[) . 

Proof of (|5.8p . We derive (|5.8p from the linearized system (|2.2|) to avoid 
working on the parameter dependent fluid domains and wave profiles. Note that 
(|2.2[) describes the evolution of the first order variations of the wave profile and 
the stream function in the travelling frame of the basic wave. The basic wave is 
(rj e (x; c), 4) e (x, y; c)J in its travelling frame (x + ct, y, t). Here, the stream func- 
tion ip e and the relative stream function tp e are related by ipe — ipe — cy, and thus 
ip e — > when |x| — > oo. As an example to illustrate the ideas, we first consider 
a perturbed solution with a trivial translation (ry e (a; + 8x; c), V> e (a; + Sx, y; c)). 
The first order variations Sx {r\ e x^ex) and SxP ex satisfy the linearized system 
(|2.2p . so we have 

=0 in T> e , 
d d - 

-r(Vey11ex) + ~rVex = On S e J 

eta; da; 

-Pex + fej/^ez =0 on 5 e ; 

-^(.WeyOnWex) + ^ 

^ e:r =0 on {y = -h} 



(ip ey d n 4> ex ) + —P ex = on S e 



from which we get A^tpex (x) = 0. Since ip ex = ipex, this recovers A°ip e x = 
which is proved in Lemma 14. II The solitary wave with the speed c + 5c is given 

by 

(r] e (x + Set; c + Sc),tp e (x + Set, y;c + 5c)) 
in the (a; + ct, y, t) frame. The first order variation are 

{r) e {x + Set; c + Sc),ip e (x + Sct,y; c + 8c)) - (r] e (x; c), ipe(x,y; c)) 
= 8c (r) ex (x;c)t + d c n e (x; c),-ip ex (x,y; c)t + d c ip e (x, y ; c)) . 

So 

(r) ex (x; c)t + d c n e (x; c), ip ex {x, y; c)t + d c ip e (x, y; c)) 

and P ex (x, y; c)t + d c P e (x, y; c) satisfy the linearized system (12.21) . By using the 
above linear system for (r) ex ,ip ex ), we get 

Adc'ipe = in V e , 

Vex + ^(^eydcVe) + ^ (d c i; e ) = on S e ; (6.2) 
d c P e + P ey d c r] e = on S e ; (6.3) 

dn^ex + ^(^eyd n (d c $ e )) + -^j- (8 c P e ) = On S e ; (6.4) 

dc'ipe = on {y = -h}. 
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From 

dcVe = —. — (d c 4>e + rje) ■ (6.5) 

tyey 



Since d n ip ex = ^ (4>ex) — ^ ( u e) and u e — > when |x| — * oo, combining (|6.5 
with (|6.4|) and (|6.3|) . we have 



9 n (d c 1p e ) = ^— (u e - PeydcTJe) = 7^ ~ -7^ ~ ^-d c Tp e . (6.6) 

So 



P _ 

Vey 



'ey 


Yey 


U e 


Pey^le 


Ipey 


1b 2 

Vey 











that is, 

4>ey 

Lastly, we show that at a turning point /Lto, A^d^ipe = 0. Indeed, the first order 
variation of 

(rj e (x + Set; no + 6/i), ip e (x + Set, y; fj, + Sfi)) - (r) e (x; [j, ),ip e (x, y; fi )) 

is Sfj,{d fl r] e ,d l j,'tp e ), since 5c = O (\5fj,\ 2 J is of higher order. So (d^,r] e , fl^ e ) 

satisfies the same system for (r]ex,'4>ex), which yields JPjd^e = 0. H 
Proof of (15791). We have 



Since 



P(c)= I r] e — (<p e ) dx = j r] e (x;c)d n ip e ((x,r) e (x;c);c))dx 



^ (d n 1pe (x, T] e (x] c) | c)) 

= -J- (ipey (X, T]e (x\ c) J c) — T] ex {x\ c) Ipex (x, T) e (x\ c) | c)) 
= d y (d c %p e ) + lpeyyd c r] e - 8 X (d c 7] e ) 1p 

ex >fex (d x (9 c 1pe) + lpexyd c r] e ) 

= (d y (dc'tpe) - r\ ex d x (d c i>e)) - d c rj e (ip exx + r) ex ip eX y) - d x (d c r) e ) ip ex 
= d n (d c tp e ) - dcVe-^; (ipex) - (dcVe) 4>ex = d n (d c $ e ) - — (d c rj e ip ex ) 
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we have 



j<9 c 77e<9„V> e + 7? e (d n (d c -ip e ) - (d c n e ipex)^ I da; 
{{d c r] e (u e - i] ex ip ex ) + n e d n (d c 4> e ) + r]exd c T] e ip ex } dx 
{u e d c r) e + T] e d n (d c ip e )}dx 

fey I \ fey fey fey 

(by equations (|6.5p and (|6.6|) ~) 



^> 5^ + ~^Ve ) ~ (dcje, + -f^ £ 

fey fey I \ fey fey 
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